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Abstract 

The purpose of this note is to show that the cross section factorization relation 
o~nn{s) I V 7 p(s) = a^ p {s) / a 11 {s) is satisfied experimentally in the energy domain 8 < 
y/s < 2000 GeV, where the cr's are total cross sections and nn denotes the even portion 
of the pp and pp total cross section. A convenient phenomenological paramaterization 
for a global simultaneous fit to the pp, pp, jp and 77 total cross section data together 
with the p- value data for pp and pp is provided by using real analytic amplitudes. 
Within experimental errors, we show that factorization is satisfied when we unfold the 
published 77 data which had averaged the cross sections obtained by using the two 
different PHOJET and PYTHIA Monte Carlo results. Our analysis clearly favors the 
PHOJET results and suggests that the additive quark model, together with vector me- 
son dominance, allows one to compute 0" 7p (s) and cr 77 (s) from o~ nn {s) with essentially 
no free parameters. The universal p-value predicted by our fit, i.e., p nn = p lp = p~ n , 
is compared to the p-value obtained by a QCD-inspired analysis of pp and pp data, 
including the p-air cross sections from cosmic rays. The p- values obtained from the two 
techniques are essentially indistinguishable in the energy region 8 < y/s < 2000 GeV, 
giving us increased confidence in our parameterization of the cross sections needed for 
the factorization relation. 
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In this note we investigate experimentally the cross section factorization relation 



(Tjp(s) ^(s)' 

where the ex's are the total cross sections and a nn , the total nucleon-nucleon cross section, is 
the even (under crossing symmetry) cross section for pp and pp scattering. 

Using eikonals for 77, ^p and the even portion of nucleon-nucleon scattering, Block and 
Kaidalov[l j have proved the factorization relation of eq. (0) by assuming that the ratio of 
elastic scattering to total scattering is process-independent, i.e., 



( 0" e l as ti c (sj \ / (^elastic (s) \ / Cr o i as ^j c (s 



V o to t(s) / V <7 tot (s) J \ a tot (s 



77 \ »""\ / / 7 p 

They have further shown that 




for all s. (2) 



Pnn(s) = P 7P (s) = p 77 (s), (3) 

where p is the ratio of the real to the imaginary portion of the forward scattering ampli- 
tude. These theorems are exact, for all s (where \fs is the c.m.s. energy), and survive 
exponentiation of the eikonal (see ref. £Q). The assumption of eq. (J2J) implies that all of 
the processes approach a black disk in the same way. In the Regge approach, factorization 
breaks down in general, for singularities other than a simple pole in the complex angular 
momentum plane. However, since the radius of interaction for the Pomeron exchange which 
determines an eikonal grows with energy as R 2 = R 2 + a' P ln(s/so) and the Pomeron slope 
a' P is very small, it was argued by Block and Kaidalovfl] that the factorization relations of 
eq. (1) were expected to be valid to a good accuracy, even in a Regge model. Indeed, one 
can give "factorization-like" relations for the residue-like constants associated with double 
or triple poles in the angular momentum plane from t-channel unitarity with certain extra 
provisos [see Cudell et. a/.[S]]. 

The strategy of this paper is to test factorization ( eq. (JTJ)) empirically by making a global 
fit to all of the experimental data for pp, pp, jp, and 77 total cross sections and the pp and 
pp p- values, i.e., making a simultaneous fit to all of the available experimental data using the 
factorization hypothesis (along with a minimum number of parameters), and seeing if the 
X 2 to this global fit gives a satisfactory value. We find that a convenient phenomenological 
framework for doing this numerical calculation is to parameterize the data using real analytic 
amplitudes that give an asymptotic In 2 s rise for the total cross sections[2j, for which there 
are mounting evidences for phenomenological success compared to other forms such as a 
power form (see for example the COMPETE collaboration [3] ). We then make the cross 
sections satisfy factorization and test the value of the overall x 2 to see if the factorization 
hypothesis is satisfied. 

We will show that the factorization relation (Tnn ^ s \ = ^mi£i j s satisfied experimentally 
when we use the PHOJET Monte Carlo analysis of the 77 cross section data, rather than 
the published values [13 Ej. We emphasize that the fit using real analytic amplitudes is only 
phenomenological and is used to provide a convenient analytical structure for the comparison 
of the shapes of the total cross sections for the three processes as a function of the energy. 
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The COMPETE collaboration J3J has also done an analysis of these data, using real 
analytical amplitudes. However, there are major differences between our analysis and the 
one done by the COMPETE group. In order to test factorization, 

• we fit simultaneously pp, pp, jp and 77 data assuming complete factorization using 
the same shape parameters, whereas they fit each reaction separately, using different 
shape parameters 

• we fit individually the two cr 77 sets of L3^3] and OPAL^3] data that are obtained 
using the PHOJET and PYTHIA Monte Carlos and do not use their average (the 
published value quoted in the Particle Data Group compilations), since the two sets 
taken individually have very different shapes and normalizations compared to their 
experimental errors. We emphasize that this individual fitting of the 77 data, i.e., a 
detailed understanding of the experimental situation, is key to our analysis. 

At the end of our computation, we investigate whether the overall \ 2 is satisfactory. 

Using real analytic amplitudes, we calculate the total cross sections cr nn , o~ lp and a 77 , 
along with the corresponding p-values. This (numerically convenient) technique has a hal- 
lowed tradition, being first proposed by Bourrely and Fischer |1| and utilized extensively by 
Kang and Nicolescujn] and more recently by Block and collaborators [HI Hj by Block[5], by 
Block and Pancheri|10). and by the COMPETE group 3J. This work follows the procedures 
and conventions used by Block and Cahnjnj. The variable s is the square of the c.m.s. en- 
ergy, p is the laboratory momentum and E is the laboratory energy. We will use a scattering 
amplitude that gives a total cross section that rises asymptotically as ln 2 (s). In terms of the 
even and odd forward scattering amplitudes / + and /~ (even and odd under the interchange 
of E — > — E), the even and odd total cross sections a even and 0" o dd are given by the optical 
theorem as 

o + = — 1m/ and cr_ = — Imj . (4) 
p p 

Thus, 

o~p P = °~+ + V- and a pp = a + — o_ (5) 
The cross section a nn , referred to in the factorization theorem of eq. (0), is given by 
An 

o- nn = — Im/ + , (6) 
p 

i.e., the even cross section. The unpolarized total cross sections for jp and 77 scattering 
are, in turn, given by 

47T 4:71 

In all of the above, the /'s and cr's are functions of s. 
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We further assume that our amplitudes are real analytic functions with a simple cut 
structure [U]. We will work in the high energy region, far above any cuts, (see ref.0, p. 587, 
eq. (5.5a), with a = 0), where the amplitudes simplify considerably and are given by 



47T 



and 



—f + (s) =i{A + (3[ln{s/s ) - iti/2] 2 + cs^e^ 1 ^ 2 } 



P 



(8) 



(9) 



where A, /3, c, sq, D, fi and a are real constants. We ignore any real subtraction constants. 
In eq. (JHJ), we have assumed that the nucleon-nucleon cross section rises asymptotically as 
In 2 s. Using equations along with eq. (jEJ) and eq. (jHJ), the total cross sections a pp , a pp 
and o nn for high energy scattering are given by 



<7 Pp (s) = A + (3 



In 2 s/sq — 



TV 



+ c sin(7r y u/2)s^ 1 - L>cos(7ra/2)s a " 1 , 



(10) 



a pp (s) = A + (3 



In 2 s/sq — 



7T 



+ c sin(7r y u/2)s^ 1 + D cos(7ra/2)s a_1 , 



(11) 



a nn (s) =A + p 



Ti- 



ll! S/S - — 



+ c sin(7r/i/2)s 



H-l 



(12) 



and the p's, the ratio of the real to the imaginary portions of the forward scattering ampli- 
tudes, are given by 



Pp P {s) 



p PP {s) 



Pnn{s) 



f3iT\ns/so - c cos(vr/i/2)s^ 1 - Dsm(7ra/2)s 0-1 



a 



(3 7T In s/s - c cos(tt / u/2)s^ 1 + D sm('Ka/2)s a ~ 1 



a 



/37rlns/s — c cos(7r///2)s 



CTr, 



(13) 
(14) 
(15) 



If we assume that the term in c is a Regge descending term, then \x = 1/2. 

To test the factorization theorem of eq. ((TJ), we write the (even) amplitudes f 1P and / 77 

as 



and 



Air 
V 



An 



/ 7p (s) = iN{A + /3pn(s/s ) - m/2} 2 + cs^ 1 ^ 1 -^ 2 } , 



—f^(s) = iN 2 {A + P[\n(s/s ) - in/2] 2 + cs^e^ 1 '^ 2 } 



(16) 



(17) 
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where N is the proportionality constant in the factorization relation ann ) s > = f^sifi — jv. 

1 1 ■ a-ypys) a"fy(s) 



We note, from eq. (jEJ), eq. (fT6|) and eq. (fTTjl . that 

/5 7r In s/sq — c cos(7r/i/2)s M_1 



A + /3 (In 2 s/s - f) + c sin(7r/V2)s"- 



automatically satisfying the Block and KaidalovpQ relation of eq. Q. 

In the additive quark model, using vector dominance, the proportionality constant N = 
fi^had' wnere -Phad * s t ne probability that a photon turns into a vector hadron. Using (see 

Table XXXV, p.393 of Ref. HQ) § = 2.2, § = 23.6 and § = 18.4, we find 

« ^v^r = V249, (19) 

where = p, tu, <fi. In this estimate, we have neither taken into account the continuum vector 
channels nor the running of the electromagnetic coupling constant, effects that will tend to 
increase P^ ad by several percent as well as give it a very slow energy dependence, increasing 
as we go to higher energies. In the spirit of the additive quark model and vector dominance, 
we can now write, using N = in eq. (fTB]) and eq. (JT7J), 



In s/s - — 



+ c sin(7r/i/2)s^ 1 (20) 



and 



In 2 s/sq 



7T 2 



+ c sin(7r / u/2)s^ 1 (21) 



with the real constants A, (3, s , c, D and P^ad being fitted by experiment (assuming a = 
H = 1/2). One might choose to vary the Regge intercepts \x and a in the fits. Since we 
want to test the goodness of the cross section factorization relations, we need a reasonable 
rendition of the even hadronic amplitude, i.e., a plays no role. Further, a small deviation of 
\i from 1/2 also gives no significant difference to our fit to shape of the hadronic data in the 
energy region of our interest, which is explicitly supported by the results of the COMPETE 
Collaboration [3J. Total cross sections for 7p scattering have been measured for energies up 
to ~ 200 GeV, while total cross sections for 77 scattering from the OPAL^] and the L3[13 
collaborations are now available for c.m.s. energies up to ^130 GeV. 

In fitting the 77 data, one might be tempted to use the 77 cross sections — along with 
their quoted errors — that are given in the Particle Data Group [12 cross section summary. 
However, on closer inspection of the original papers, it turns out that results quoted by the 
PDG are the averages of two independent analyses performed by both the OPAL ^3] and 
L3[T3*] groups, using the two different Monte Carlo programs, PHOJET and PYTHIA. The 
error quoted by the Particle Data Group was essentially half the difference between these two 
very different values, rather than the smaller errors associated with each individual analysis. 

The Monte Carlo simulations used by OPAL and L3 play a critical role in unfolding 
the 77 cross sections from the raw data. To quote the OPAL authors [Hj. "In most of 
the distributions, both Monte Carlo models describe the data equally well and there is no 
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reason for preferring one model over the other for the unfolding of the data. We therefore 
average the results of the unfolding. The difference between this cross section and the results 
obtained by using PYTHIA or PHOJET alone are taken as the systematic error due to the 
Monte Carlo model dependence of the unfolding." For the testing of factorization, there is 
good reason for possibly preferring one model over another, since the two models give both 
different normalizations and shapes, which are vital to our analysis. Hence, we have gone 
back to the original papers f3J E] and have deconvoluted the data, according to whether 
PHOJET or PYTHIA was used. These results are given in Fig. ^ Clearly, there are major 
differences in shape and normalization that are due to the different Monte Carlos, with the 
PYTHIA results significantly higher and rising much faster for energies above ~ 15 GeV. 
On the other hand, the OPAL and L3 data agree within errors, for each of the two Monte 
Carlos, and seem to be quite consistent with each other, as seen in Fig. [TJ 

For these reasons, we will make three different fits, whose results are shown in Table ^ 
Fit 1 is a simultaneous x 2 fit °f ec L- (HHJ), (fTTj) . (fT3")) . ffT^j) and (}2T)j) to the experimental a pp , o pp , 
Ppp, Ppp and a w data in the c.m.s. energy interval 10 GeV < yfs < 1800 GeV, i.e., we don't 
include the 77 data. We next make two different simultaneous x 2 fits of eq. (fTUj) . (fTTj) . (fT3"j) . 
(fTH) . (|2U|) and (J2T)) to the experimental a pp , a pp , p pp , p pp , a lp and the unfolded <7 77 , using 
either PHOJET or PYTHIA results, in the c.m.s. energy interval 10 GeV < < 1800 
GeV. Fit 2 uses a 77 from PHOJET unfolding and Fit 3 uses <r 77 from PYTHIA unfolding. 
In order to account for possible systematic overall-normalization factors in the experimental 
data, the cross sections for L3 are multiplied by the over all-renormalizat ion factor Al3 and 
those for OPAL are multiplied by an overall-renormalization factor Aopal, with these factors 
also being fitted in Fits 2 and 3. 

^From Fits 1, 2 and 3, we see that the major fit parameters A, [3, so, D, c and P^ d are 
the same, within errors. The purpose of Fit 1 was to show the robustness of our procedure, 
independent of the 77 data. 

However, when we introduce the unfolded 77 cross sections in Fits 2 and 3, we see that the 
results strongly favor the PHOJET data of Fit 2. The x 2 /d.f. jumps from 1.49 to 1.87 (the 
total x 2 changes from 115.9 to 146.0 for the same number of degrees of freedom). Further, 
and perhaps more compelling, the normalizations for both OPAL and L3 are in complete 
agreement, being 0.929 ± 0.037 and 0.929 ± 0.025, respectively. Thus, they differ from unity 
by ~ 7 ±3%, compatible with the experimental systematic normalization error of 5% quoted 
by L3 . The PYTHIA results from Fit 3 have normalizations that disagree by ~ 14% and 
~ 19% for OPAL and L3, respectively, in sharp disagreement with the 5% estimate. Thus, 
from here on, we only utilize the PHOJET results of Fit 2, whose parameters are given in 
Table 1. 

Using the parameters of Fit 2, we find that P^ad = 1/(233.1 ± 0.63), which is in reason- 
able agreement with our preliminary estimate of 1/249, being 6% larger, an effect easily 
accounted for by continuum vector channels in 7^ reactions that are not accounted for in 
the estimate of eq. (fT9"j) . 

The fitted total cross sections a pp and a pp from eq. (fTUj) and eq. (fTTj) are shown in Fig. 
I2J along with the experimental data. The fitted p-values, p pp and p pp from eq. (|T3"|) and 
eq. ()14|) are shown in Fig. along with the experimental data. The fitted total cross 
section a 7P = |-Ph ad o" nri from eq. (J2T?)) is compared to the experimental data in Fig. HI using 
-^had = 1/233. The overall agreement of the pp, pp and 7^ data with the fitted curves is 
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quite satisfactory. We now turn our attention to the 77 data. 

The fitted total cross section <t 77 — (l-^had)^^™™ from eq. (|21j) is compared to the exper- 
imental data in Fig. again using P^ ad = 1/233. The experimental data plotted in Fig. 
are not renormalized, but are the results of unfolding the original experimental results, i.e., 
use iVopAL = A^L3 = 1. We see from Fig. El that within errors, both the shape and normal- 
ization of the PHOJET cross sections from both OPAL and L3 are in reasonable agreement 
with the factorization theorem of eq. (|I|L whereas the PYTHIA cross sections are in distinct 
disagreement. This conclusion is born out by the x 2 's of Fit 2 and Fit 3 in Table 1. 

Finally, the fitted results for <t 77 , using the parameters of Fit 2, are compared to the 
renormalized OPAL and L3 (PHOJET only) data in Fig. |U1 The agreement in shape and 
magnitude is quite satisfactory, indicating strong experimental support for factorization. 

For completeness, we show in Fig. Qthe expected p- value for the even amplitude, from 
eq. (JTHJ). Also shown in this graph is the predicted value for p nn found from a Q CD-inspired 
eikonal fit by Block et al. ^3] to pp and pp total cross sections and p- values from accelerators 
plus p-air cross sections from cosmic rays. The agreement between these two independent 
analyses, using very different approaches, with one using real analytic amplitudes with a In s 2 
behavior and the other using a QCD-inspired eikonal model in impact parameter space, giv- 
ing rise to a cross section also eventually rising as Ins 2 , is most striking. In both cases, these 
two analyses give p nn = p lp = p 77 , another factorization theorem of Block and KaidalovpQ. 

We conclude that the cross section factorization hypothesis of [1 , = , is 



satisfied for nn, ■yp and 77 scattering, if one uses the PHOJET Monte Carlo program to 
analyze <r 77 . Further, we find that the experimental data also satisfy the additive quark 
model using vector meson dominance, since 

9 

p7 

Cr ' rp _ 3 had(7nn 

^77 = (fj£d) *nn, (22) 

with k = 2/3 and P h 7 ad = 1/233. 

The assumption of Block and Kaidalov[T] in eq. (J2J) that o" e i astic (s)/o" to t(s) is process- 
independent yields another factorization theoremjTj 

B lp (s) ~ B^aY { 1 

where the 5's are the nuclear slopes for elastic scattering (the logarithmic derivatives of the 
elastic scattering cross sections <icr elastic /(it, where t is the squared 4-momentum transfer). 
For 7p processes, using vector dominance, the £>'s are the slopes of the 'elastic' scattering 
reactions 

7 + p -> V + p, (24) 

where the vector meson V is either a p, uj or <fi meson. Using the additive quark model, 
eq. (j211) implies that 

2 
3' 



B Jp (s) = K,B nn (s), where k = —. (25) 
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It has been shown by Block, Halzen and Pancheri|lb'| that a x 2 fit to the available 7p data 
gives 

k = 0.661 ±0.008, (26) 

in excellent agreement with the 2/3 value predicted by the additive quark model, again 
justifying the use of 2/3 in our fits. 

We conclude that if we determine <7 nn (s), p nn and B nn from experimental pp and pp data 
for y/s > 8 GeV, we can then predict rather accurately 

in essentially a parameter-free way, by using factorization and the additive quark model with 
vector dominance. Clearly, this conclusion would be greatly strengthened by precision cross 
section measurements of both 7p and 77 reactions at high energies. 
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o-tot ~ ln 2 (s/s ) 


Parameters 


r ii 1. 


rll z. 


rll o. 




no o" 77 


_ fv^rr, T3T1TO TT7"T 

(T 77 trom rrlUJiljl 


<7 77 trom r Y 1 rllA 


A ( I 

A {mu) 


6 1 .A ± U.ol 


d ( .1± U.o 


d l.o± u. I I 


P (JXIQ ) 


0.304 ±0.023 


0.302 ±0.024 


0.307 ±0.022 


so ((GeV) 2 ) 


34.3 ± 14 


32.6 ± 16 


35.1 ± 14 


D (mb(GeV) 2(1 ~ a) ) 


-35.1 ±0.83 


-35.1 ±0.85 


-35.4 ±0.84 


a 


0.5 


0.5 


0.5 


c (mb(GeV) 2(1 ~ M) ) 


55.0 ±7.5 


55.9 ±8.1 


54.6 ± 7.3 




0.5 


0.5 


0.5 


pi 

Miad 


1/(233.1 ±0.63) 


1/(233.1 ±0.63) 


1/(233.0 ±0.63) 


-^OPAL 




0.929 ±0.037 


0.861 ±0.050 


iV L3 




0.929 ±0.025 


0.808 ±0.020 


X 2 /d.f. 


1.62 


1.49 


1.87 


d.f. 


68 


78 


78 


total x 2 


110.5 


115.9 


146.0 



Table 1: Fit 1 is the result of a fit to total cross sections and p- values for pp and pp, along with 
a™. Fit 2 and Fit 3 are the results of fitting total cross sections and p- values for pp, pp and ct 7P , 
as well as including the <7 77 data from the OPAL and L3 collaborations. Fit 2 uses the results of 
unfolding <t 77 with the PHO JET Monte Carlo, whereas Fit 3 uses the results of unfolding <t 77 with 
the PYTHIA Monte Carlo. The overall-renormalization factors iVopAL and A^L3 are also fitted in 
both Fit 2 and Fit 3. The fitted parameters are the ones that have statistical errors indicated. 
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• L3, using PHOJET 

■ OPAL, using PHOJET 

O L3, using PYTHIA 

□ OPAL, using PYTHIA 



M si 



t/s, in GeV 



Figure 1: OPAL and L3 total cross sections for 77 scattering, in nb vs. s/s, the c.m.s. energy, in GeV. 
The data have been unfolded according to the Monte Carlo used. The solid circles are the L3 data, unfolded 
using PHOJET and the open circles are the L3 data, unfolded using PYTHIA. The solid squares are the 
OPAL data, unfolded using PHOJET and the open squares are the OPAL data unfolded using PYTHIA. 




Figure 2: The dotted curve is a pp , the predicted total cross section for pp reactions (from Fit 2), in mb, 
and the solid curve is ap p , the predicted cross section for pp reactions (from Fit 2), in mb vs. ^/s, the c.m.s. 
energy, in GeV. The circles are the experimental data for pp reactions and the squares are the experimental 
pp data. 



- 9 - 



0.3 



0.2 



a 0.0 



-0.2 



-0.3 




• PP 

■ pbar-p 

pp, log 2 (s) analytic amplitude fit 
pbar-p, log 2 (s) analytic amplitude fit 
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10000 
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Figure 3: The dotted curve is p pp , the predicted ratio of the real to imaginary part of the forward scattering 
amplitude for pp reactions and the solid curve is pp p , the predicted ratio for pp reactions, vs. *Js, the c.m.s. 
energy, in GeV, from Fit 2. The circles are the experimental data for pp reactions and the squares are the 
experimental pp data. 
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300 
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1000 



Figure 4: The curve is a lv — t^had * 111 ' ^ ne P re dicted total cross section for 775 reactions, in fib vs. \/s, 
the c.m.s. energy, in GeV, from Fit 2. The circles are the experimental data. 
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■ OPAL, using PHOJET, uncorrected 

□ OPAL, using PYTHIA, uncorrected 

O L3, using PYTHIA, uncorrected 
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10 



100 
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Figure 5: The curve is — (§-f > had) 2 "' nn ' predicted total cross section from Fit 2 for 'y^y reactions, in 
nb, vs. t/s, the c.m.s. energy, in GeV. The open squares and circles are the experimental total cross sections 
for OPAL and L3, respectively, unfolded using the PYTHIA Monte Carlo. The solid squares and circles are 
the experimental total cross sections for OPAL and L3, respectively, unfolded using the PHOJET Monte 
Carlo. 




Figure 6: The curve is avy-y — (l^had) 2 ^ 11 ' predicted total cross section from Fit 2 for "yy reactions, 
in nb, vs. y/s, the c.m.s. energy, in GeV. The squares and circles are the total cross sections for OPAL and 
L3, respectively, unfolded using PHOJET, after they have been renormalized by the factors A^opal = 0.929 
and A L3 = 0.929 found in Fit 2 of Table 1. 
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Fi gure 7: The solid curve is p 77 = p rp = p n n, the predicted ratio (from Fit 2) of the real to imaginary 
part of the forward scattering amplitude for the even amplitude vs. y/s, the c.m.s. energy, in GeV. The 
dotted curve, shown for comparison, is p nn , the result of a QCD-inspired cikonal fit|15| to pp and pp data 
that included cosmic ray p-air data. 
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